Preprint UTPT-99-01 



Infrared behavior of graviton-graviton scattering 

John F. Donoghue 1 and Tibor Torma 2 

1 Department of Physics and Astronomy, University of Massachusetts 
Amherst, MA 01003, USA* 
q^" ■ 2 Department of Physics, University of Toronto 



in 
o 

On 
On 

23 



Toronto, Ontario, Canada M5S1A7^ 



c : February 1, 2008 

d 
i— » 

oo 

<N . 

T ~~ l 1 Abstract 
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\ The quantum effective theory of general relativity, independent of the even- 

tual full theory at high energy, expresses graviton-graviton scattering at one 
loop order O (E 4 ) with only one parameter, Newton's constant. Dunbar and 
Norridge have calculated the one loop amplitude using string based techniques. 
We complete the calculation by showing that the divergence which remains 
in their result comes from the infrared sector and that the cross section is finite 
and model independent when the usual bremsstrahlung diagrams are included. 

! 1 Introduction 

> ■ 

•rH , 

The simplest low energy process in quantum gravity is graviton-graviton scattering. 
Although experimentally unobservable, this reaction forms an interesting theoretical 
laboratory that illustrates the workings of quantum gravity. If general relativity is 
the correct low energy classical theory of gravity, then its quantum theory forms an 
effective field theory capable of analyzing the low energy quantum effects. Graviton- 
graviton scattering is particularly useful in illustrating the logic of predictions in a 
quantum effective theory. Indeed, at one loop order this reaction provides a model- 
independent quantum prediction of general relativity. 
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At tree level, the graviton-graviton scattering amplitude is simple in the helicity 
basis, although the calculation to obtain this result from the Einstein action is not so 
simple. With + (— ) representing helicity +2 (— 2)[|, all tree amplitudes for 1+2 — > 3+4 
vanish except those related to A tree (++; ++) by crossing and |J 

A tree (++:++" - 



' ' 4 tu 

Here k 2 = 32nG and s = (pi + P2) 2 , t = (pi — P3) 2 , u = (pi — pi) 2 denote the usual 
Madelstam variables. 

It is simple to show that graviton-graviton scattering should be finite and pa- 
rameter independent at one-loop order ||. In the effective low energy theory [§ 
gravitational effects are expanded in a derivative expansion with all terms satisfying 
general covariance 

r r 2 1 

Sgrav = / d A X y/g — R + C\R +C 2 R^R^ + ... + £ matter (2) 

Here k 2 = 32nG and G is Newton's constant, c+2 are unknown dimensionless param- 
eters which contain information about the (presently unknown) ultimate high energy 
theory. A third covariant of order R 2 , R^vupR^ 01 ^ can be removed in four dimensional 
space-time through the use of the Bianchi identities. Since the curvature involves two 
derivatives of the metric, the Einstein action (the term with R) is seen to be of order 
E 2 , while R 2 and R fJiU R' J ' u are of order E A . 

Loop diagrams obey a power-counting theorem |J. One loop diagrams formed 
from vertices given by the Einstein action yield effects at order E A - any process 
with more loops is higher order in the energy expansion. The ultraviolet divergences 
at one loop necessarily have the same structure as the local Lagrangian in Eqn. (0), 
which means that they must be proportional to R 2 or R^R^ . Then, at this order, the 
ultraviolet divergences can be absorbed into renormalized values of the the parameters 
c+2- These renormalized constants are unknown and will be different depending on 
the nature of the theory that forms the ultimate correct high energy theory which 
includes gravity. In this sense these parameters are model dependent. However, they 
do not contribute to the process of graviton-graviton scattering. At the order that we 
are working, the R 2 Lagrangians are applied to form vertices for on-shell amplitudes, 



1 Note that in our notation crossing also requires one to flip the ± sign for the affected gravitons. 
This implies in particular that A{ ; ++) must be a symmetric function of s, t and u. 
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which is to say that the equations of motion are satisfied for the external states. 
However, the equations of motion for the purely gravitational sector are = 0, 
and hence R = also. Thus the effects of both of the R 2 terms in Eqn. (|2|) vanish in 
purely gravitational processes. It is this argument that tells us that graviton-graviton 
scattering is finite and independent of any unknown parameters at one loop order. 

The power counting theorem is manifest in the one-loop results calculated by 
Dunbar and Norridge J/J. The one-loop amplitude is formed by using the lowest 
order tree amplitude twice in order to produce a loop diagram, and hence carries 
coupling constants k 4 ~ G 2 N . Dimensionally this requires that the result carry four 
powers of the external energies. This is seen in the results: 
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and all logarithms with negative arguments are understood to have a —in imaginary 
part. Note that this represents a real tour- de- force. Done in conventional field theory, 
the calculation is formidably difficult. It is a tribute to the string based techniques 
that the results are obtainable with less then Herculean effort. Indeed, after cal- 
culating the graviton loops, the authors write down the result for massless scalars, 
fermions and photons in the loops in just a few lines. However, the result is not tied 
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to the validity of string theory as a fundamental theory - the technique is simply an 
efficient way to calculate the results of usual (quantum) general relativity. 

One notices that the one loop amplitude in Eqn. @ contains a factor of -, i.e. 
it is not finite. At first, this seems to contradict the general reasoning given above. 
However, in the complete calculation of the physical process of graviton scattering, 
there will also be bremsstrahlung diagrams describing the radiation of soft gravitons 
off the external graviton lines. When calculated in a d = 4 — 2e dimensional phase 
space these infrared effects also bring in a \ factor. If the divergence in Eqn. is 
an infrared divergence, and if the effective field theory of gravity behaves as a proper 
effective field theory, then the infrared loop effects should be canceled against the soft 
radiation. While there are good reasons for believing that the gravitational effective 
field theory should be well behaved in the infrared, the long-standing doubts about 
quantum gravity make it worthwhile to check this property in the only complete 
calculation available. In Ref. || it was shown that the scattering of spin-0 fields 
is infrared finite even in the limit when their masses vanish. However, it was only 
conjectured there that the same is true for massless matter of higher spin (a situation 
similar to graviton-graviton scattering). One also notes that the scale of the logarithm 
is not defined. This is an indication that the calculation is incomplete. We will see 
that the scale in the logarithm comes from an infrared regulator for soft gravitons. 
Finally, part of our motivation comes from a minor quibble with the argument given 
above. In the effective Lagrangian we removed the (R^ap) 2 term by the use of an 
identity that is only valid in exactly four dimensions. Indeed, in any higher dimension 
the argument given would not apply, and the graviton scattering amplitude would 
contain a model dependent parameter. This means that in the quantum theory we 
can only be certain of the result if we use a regularization scheme that works in 
four dimensions. However, the only scheme that we know about that preserves the 
symmetries of general relativity is dimensional regularization, and it was that scheme 
used in Ref. 0. While it is unlikely that the regularization scheme would lead to 
an extra divergence, we also want to confirm that the residual divergence is not an 
artifact of the ultraviolet regularization. 
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2 Soft gravitons in graviton-graviton scattering 

We will explicitly calculate the divergences in the one-loop differential cross section for 
graviton-graviton scattering. We will find a complete cancellation of infrared diver- 
gences when we calculate the cross section up to O(kP), including the Bremsstrahlung 
graphs, as shown in Fig. [I]. 



Figure 1: The expansion of the cross section in k in graviton-graviton scattering. The 
quantity A tree represents the sum of all tree level diagrams. Solid lines represent hard 
gravitons, wavy lines are soft gravitons. 

In this figure we explicitly show all factors of k. The first term in the figure (and 
five additional graphs, not shown, with graviton exchanges between various pairs of 
external legs) is already included in the full one-loop scattering calculation and has an 
infrared divergence. This divergence is canceled by another divergence in the second 
term in the figure, a soft Bremsstrahlung process, which should be added as it is 
degenerate in energy with pure hard scattering. The second line in the figure shows 
that the actual cancellation occurs in the 0(k, 6 ) terms because the leading 0(/t 4 ) is 
tree level and infrared finite. 

We will derive a general formula for the infrared divergences which uses the on- 
shell Born amplitude. The most convenient regularization procedure is dimensional 
regularization. We calculate the IR divergent part of the graviton radiation term in 
Fig. [I] and show [cf. Eqn. (|26|)1 that to do so we only need to know the on-shell tree 
level amplitudes^4 tree (Ai, A 2 , A 3 , A 4 ). We always work in the helicity basis and A; = ± 
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stands for the helicity of the i hard graviton. The only divergence occurs when 
the gravitons have helicity assignments ++; ++ (and in the cases related to this by 
crossing) and that is the only case when the tree level amplitude is nonvanishing. In 
the following we show [Eqn. fll3|) l that soft graviton radiation does not flip the hard 
graviton spins so that all IR divergences are proportional to the tree amplitude with 
the same helicity. 

The amplitude with one soft graviton radiation is the sum of the four diagrams 
in Fig. 0. We first calculate the contribution from 2(a) 




Figure 2: The four Feynman diagrams that contribute to soft graviton radiation at 
lowest order in hard graviton-graviton scattering. 



A$ = (1, 2; 3, h + k) - 1 — P^, v , i^ Xp ef & (5) 



where writing a number in the argument of the tree amplitude means putting those 
lines on shell and multiplying by the appropriate & v polarization tensor. The matrices 
/ and P denote 

I a p,yS = 2 {Va-yVpS + VaSVp-y) (6) 

and 

The gauge invariance of the tree level amplitudes implies, for k — > 0, 



k1A^ e (l,2;3,h + k)=0(k) (8) 

and also 

V ^A flu (l,2;3,h + k) = 0(k). (9) 
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These restrictions can be derived as follows. Gauge invariance implies that the on- 
shell amplitude is unchanged under shifting the polarization tensor by 

— > + {k"C + -k-i rT) , (10) 

with any four-vector a transformation that keeps k^e^ zero. In order for an 
amplitude A^e^ to be invariant under such replacement, we need for any on shell 
momentum k 

2UKAn = (k-0^. (11) 
This must hold for any £, hence we have Eqns. (§,0). 




i ft P 1 rv/9 



K " T a^s( k uk 2 ,h) 



Figure 3: The graviton propagator and the triple gluon vertices in harmonic gauge. 
For an expression of r^g 7<5 see Eqn. (0). 

The graviton propagator and the triple graviton couplings are shown in Fig. ^| 
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Putting together Eqns. (^| [T^) , we arrive at a simplified expression 



Arad _ Ke^(k,\)k% .trees-, o o a\ i KO(k) . . 

where the "IR" index emphasizes that we keep only the leading term when k — * 0. We 
observe indeed that the O (j^J term is proportional to the Born amplitude without 
flipping any of the hard particles' spins. 

Now we add on the contribution from Fig. [|b,c,d. The result is 



jrad __ _ K j^ Tee (\ 2 3 4) V" k^e^k, \)k^ ^ k Q(k) 

n=l V^n Vnk) 2 + Z€ k ■ k n 

Next we need to square this amplitude and sum over the soft graviton spin: 

4 k^k u T\ a(k)k a k lS 
\^ I Arad\2 _ 2 i A tree\2 i "j 1 V',a/H / V ftj ft-j , . 

V irI " 1 1 ^PWWTW' ( 5) 

where the sum over graviton polarization tensors is 

u^' al3 (k) = ]T ^±(k)^±(k)7^(k)ei(k) = - (n" a ir /J + n^rr" - n^n a/3 ) (ie) 

and 

IF"(fc) = k"X u + k v \» - (k - A) 77^ (17) 

with an arbitrary vector A, same for all terms in the sum, chosen as = (1,0). 
Doing the algebra in the above formula gives us 

Arad\2 \A tree \ 2 ^ (cos7y - cos a^ cos a,,) 2 - i sin 2 «j sin 2 ttj 

^ 4/r ^ (1 — cos ay (1 — cos ay 

where k now stands for the energy of the soft gluon (not the four- momentum), and 
7jj is the angle between the (d— l)-dimensional momenta of the hard gravitons, on is 
the angle between the i th hard and the soft gravitons; Ei is the CM energy of the i th 
graviton and rji — +1 (—1) for incoming (outgoing) hard gravitons. 

At this point we make a comment on how dimensional regularization works. In 
the one-loop amplitude we find a - divergence in dimensional regularization. As 
pure gravity is one-loop finite and all the divergences in one-loop graviton-graviton 
scattering come from the pure gravity part only, all of this - should be of infrared 
origin and consequently be canceled by the square of the amplitude A rad . However, 
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jrad j^ se jf j s a ^ ree j eve i amplitude which does not diverge; the canceling - factor 
comes from the phase space integral. One might wonder then why we are not getting 
too much divergence: the leading term is -k, so the phase space integral introduces 



which is logarithmically divergent. In the same time the angular integration is also 
divergent and we find that dimensional regularization does not handle correctly an 
integral of the type 



(here a is the angle between the direction of n and a fixed direction.) The above 
Euler function is 



finite, although the integral includes a severe collinear singularity. Fortunately, in 
our case, we will not encounter this problem: the spins "conspire" so that there is an 
additional angular factor which takes away all collinear singularities in this integral. 
In other models, however, like one with elementary massless scalars, this might be a 
problem which requires further treatment. 

Now we calculate the differential cross section in d = 4 — 2e dimensions. We 
focus on the infrared region only, integrating up to a cutoff A <C ^Js and neglecting 
momenta of order A and above. Such soft graviton radiation should (and will) be 
sufficient to cancel the IR divergences due to one-loop integrals. In particular, we do 
not consider hard collinear gravitons. The divergences due to hard collinear graviton 
radiation (i.e. when one of the ctj's is small) are not canceled by loops. However, these 
divergences are all proportional to A so can be unambiguously separated from the 
soft divergences. Some rather tedious algebra leads to an integral over the direction 
m of the soft Bremsstrahlung graviton 




(19) 




(20) 




(21) 




(22) 
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The k integral has a \ infrared divergence. All divergences that are collinear and 
infrared simultaneously should come from the second integral. However, we observe 
that the numerator in the angular integral vanishes when the denominator does, 
actually canceling out the singularity. This fact is necessary to allow us to consistently 
separate collinear divergences from soft ones. In order to find the divergent part, we 
need to calculate only the leading term in 



^ (0) ( 7 ) + e^ ij (7) + --- 



dfi d _i(m) 



(cos 7y — cos ctj cos aj) 2 — | sin 2 «j sin 2 aj 
(1 — cosai)(l — cos ay) 



Substituting this into Eqn. 

] rad 

rffi d _x(n) 



we find 
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The result in four dimensions is 
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With this, we finally find for the cross section 
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(23) 



(24) 



(25) 



(26) 



We have found that the infrared divergent part is indeed proportional to the square of 

the Born amplitude. Because in the (++; ) and (++; H — ) helicity cases there is 

no IR divergence to cancel, the vanishing of the Born terms makes sure none emerges 
in the radiative process. 

In the (++; ++) helicity case we need to use the cross section formula 



( d °l99 -> gg] 
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2Re(A tree A 1 - loop ) 
(2tt) 2 2 5 s 



(27) 



in order to calculate the 0(k 6 ) contribution to the cross section (see Fig. |T|.) Using 
the Dunbar- Norridge 1-loop amplitude amplitude, Eqn. @, we find the C(k 6 ) 
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contribution to the cross section for the 2 — > 2 process: 
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We observe that the - divergence cancels when we add together Eqns. ([26|) and (|28f) . 
The finite term in Eqn. fl28|) contains an undetermined scale due to the logarithm of 
s. The occurrence of such a scale is a common feature of dimensional regularization 
in the presence of infinities. This scale is provided by the "ultraviolet" cutoff in the 
radiative cross section. Our final result for the sum of the cross sections is 
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In this form all divergences are canceled and all logarithms are dimensionless. There 
is a logarithmic dependence on the scale where we cut off the non-infrared radiative 



gravitons. We remind the reader that the finite functions / 



and 



are respectively given in Eqn. (£|) and can be extracted from Eqn. (|23|). 



3 Conclusions 

This has been the first explicit investigation of the infrared properties of a one loop 
amplitude in quantum gravity. We have achieved our goal in demonstrating that 
the effective theory of gravitation is not plagued by infrared divergences, its soft 
divergences even cancel in the case of one-loop graviton-graviton scattering, and also 
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demonstrated that, similarly to the case of QED, summation over degenerate states 
in the final state suffices to get a final and sensible cross section. 

The result for graviton-graviton scattering to one-loop order is beautiful and sig- 
nificant because it forms a low energy theorem for quantum gravity. No matter what 
the high energy theory of gravity may turn out to be, and independent of the mas- 
sive particles in the theory, as long as the low energy limit leads to general relativity 
the scattering rate must have the model independent form shown in Eqn (|29|). As 
expected, the quantum effective field theory of general relativity is well-behaved in 
the infrared. 
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